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We study the effect of an environment consisting of noninteracting two level systems on 
Landau-Zener transitions with an interest on the performance of an adiabatic quantum 
computer. We show that if the environment is initially at zero temperature, it does not 
affect the transition probability. An excited environment, however, will always increase 
the probability of making a transition out of the ground state. For the case of equal 
intermediate gaps, we find an analytical upper bound for the transition probability in 
the limit of large number of environmental spins. We show that such an environment 
will only suppress the probability of success for adiabatic quantum computation by at 
most a factor close to 1/2. 
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1. Introduction 

Theories of open quantum systems have gained renewed attention in recent years, 
thanks to their important role in quantum computation. Such theories usually study 
quantum evolution of a system in contact with an environment, which in general can 
be made of bosons^, fermionJ21 ; or localized spinJ^. At the end of the evolution, the 
properties of the system are extracted by averaging over the environmental degrees 
of freedom. 

Significant progress has been made in understanding the effect of environment on 
quantum coherence and quantum computatiorF^^. There is evidence that certain 
quantum computational models are less sensitive to loss of ph ase co herence than 
others. One of those, is adiabatic quantum computation 

(AQcj^HH. i n AQC, the 

system starts in a known ground state of a Hamiltonian Hi . The Hamiltonian then 
slowly evolves to a final Hamiltonian Hf, whose ground state encodes the solution 
to the problem. In the absence of an environment, if the evolution is slow enough, 
the system will end up in the ground state of Hf with probability close to one. The 
total time dependent system Hamiltonian is usually written as a linear interpolation 

Hs(t) = [l-s(t)]Hi+s(t)H f , (1) 

'Currently at UBC mathematics department, 
t Currently at MIT physics department. 
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where s(t) S [0, 1] is a monotonic function of time t. 

In the absence of an environment, the main source of error for AQC is non- 
adiabatic transition out of the ground state. In the small gap regime, it is closely 
related to the Landau-Zener (LZ) transitions^^, which occur when two energy levels 
pass an avoided crossing (anticrossing) . In the original LZ problem, a 2-level system 
evolves via a Hamiltonian H— — (At x +vtT z ) /2, where t x > z are the Pauli matrices. 
The spectrum of this system has an anticrossing at t—0 with gap A. If at t— >— oo 
the system is in its ground st ate, the probability of finding it in its excited state at 
time t— >oo is exactly given by^ 

P LZ =e-* A ' 2 ^. (2) 

In AQC, the system has many energy levels, and the running time is finite; 
therefore the above problem does not directly apply. However, when the gap g m 
between the first two levels is much smaller than other relevant energy scales, the 
transition probability is given by ^ to a very good approximation. In such a case, 
v oc l/tf, where tf is the running time. In order to have small excitation probability, 
one needs a long tf (oc 1/g,^). Therefore, problems with small gap are the hardest 
to solve using AQC. 

Recently, there has been compelling evidence that localized two level sys- 
tems p lay an imp ortant role in removing phase coherence in solid state quantum 
systems^EEE^nSl Such 2- level sy stems can be nuclear spind^l, magnetic impurities, 
or other 2-level fluctuator J^E^l The effect of a bosonic environment on LZ tran- 
sition has been studied extensively^SIISnS and some studies have been done for a 
spin environment!^^. In this article, we study a spin environment focusing only on 
small gap problems, relevant for AQC. 



2. Hamiltonian 

As usual we write the total Hamiltonian as H = Hs + Hb + -Hint , which includes 
system, bath, and interaction Hamiltonians respectively. If the minimum gap g m is 
much smaller than the separation of the two crossing levels from the other energy 
levels, then the slow evolution of the system close to the anticrossing will be re- 
stricted only to those two levels. Using a new coordinate e = 2E(s—s*), where E is 
an energy scale characterizing the anticrossing and s* is its position, one can write 
a two-state Hamiltonian: 

H s = -(eT z +g m r x )/2, (3) 
and the gap between the first two states is well approximated by 

.9 = Vt 2 + g 2 m - (4) 

with t XjZ being the Pauli matrices in the two-state subspace of the central system. 
We only focus on linear time sweep for which e = vt with v = 2E/tf. In this case, 
the two-state problem becomes a LZ problem. 
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We incorporate the environment via the Hamiltonians 



Hp 



1 m 



3=1 



1 m 



(5) 



Here, a x - z are Pauli matrices for the environmental spins, Bj=[BJ, Bj , J3J] is a local 
field applied to the j'-th environmental spin, Jj is the coupling coefficient between 
j-th spin and the system, and m is the number of environmental spins. Notice 
that we only take into account longitudinal coupling between the effective two- 
state system (r z ) and the environmental spins (cr?). As is shown in the appendix, 
in the case of adiabatic Grover search problem a general coupling between 
individual qubits within the system and the environmental spins leads to only a 
longitudinal coupling between the effective two-state system at the anticrossing and 
the environmental spins in the limit of large number of qubits. The same behavior is 
however expected for more general Hamiltonians because transverse coupling to the 
environment will cause large relaxation to the ground state, enabling us to solve the 
problem efficiently by relaying only on classical relaxation (i.e., classical annealing) 
which is not expected. 

In general, the effect of the environment is to split (broaden) each energy level 
into M=2 m levels, thereby splitting the original anticrossing into M 2 avoided cross- 
ings, as shown in Fig. [T] Since we are interested in small g m , we use a perturbative 
approach: H=Hq+H' , with the perturbation Hamiltonian H'= — g m T x /2 and the 
unperturbed Hamiltonian Hq which constitutes all the terms in H except H' . To 
diagonalize Hq, we introduce a unitary transformation U by 



U = 



3=1 

a ■, 



-)I+(cos0t 



")I+ (sin6n 



ef = arccos {Slf ± B* + J 3 ) / 2^ 

nf = J(Bj* 



(fl¥)a + (B*±J 3 .) a . 



sm9 j )t z , 



After the transformation, we find 

UH U^ = E ZSl ... Sm \zsi...s m )(zsi...s m \, 



UH'U^ 



1 



n 

3 = 1 



{B^ + {BW 



E- 



1 1 - 

S1 ... sm = -r {e(-iy + - ^(-i) s ' m+ sij) +(-ir (o+ + nj)]} 



(6) 
(7) 
(8) 



3 = 1 



I 



I 
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where z and Sj are eigenvalues of (I — r z )/2 and (I — cr |)/2 respectively, and 0j = 

Let \z, k) denote the state of the system plus the environment when the envi- 
ronment is in its fc-th energy level. We use degenerate perturbation theory to find 
the gap gkk' at the point where |0, k) and |1, A') cross (see Fig. [J). Assuming that 
Hq has no other degeneracies except at the crossings, we find 

m 

9kk' = 9m Y]_ \cosdj5 SjS r. + sin6»j(l - S SjS r.) , (9) 
i=i 

where 5 SjS r is the Kronecker delta function. Now for any fc, one can easily show 
that 



M m 

E sL = ul li E H 2 + sin2 jO- - VP 



fc' = l 3=1 s>. 

■2 



(10) 



= 9n 

This is a very important sum rule which will be used in our analysis repeatedly. 



3. Environment at zero temperature 

While the intermediate gaps (JHJ) are proportional to g m , their separation is indepen- 
dent of g m . Therefore, in the limit g m — > 0, there is always a regime in which the 
intermediate gaps can be considered separated enough to justify using LZ formula, 
Pkk' = e -71 " 9 "' l 2v ', for the probability of transition at the kk'-th anticrossing. One 




Fig. 1. Energy spectrum of the system plus the environment. The arrows show paths correspond- 
ing to excitation of the system. The green line shows a path in the fast passage limit. The red and 
blue lines show two interfering paths. 
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immediately finds 

\{p k y=^A-Y d ^9kwl^\=^ a2m/2v = ^, (11) 

k' = l Ik' ) 

where Plz is given by ([2|). Equation (jllj) has an important consequence when the 
environment is initially in its ground state (k=l). The probability of ending up 
in the final excited state is the probability that all intermediate LZ excitations 
occur (straight path in Fig. QJ: P e f — YikLi Pik=Pui- Hence, a zero temperature 
spin environment does not affect LZ transition probability. Similar results have a lso 
been reported for an environment made of uncoupled harmonic osci 



4. Environment at finite temperature 

We now discuss the effect of an excited (equilibrium or nonequilibrium) environ- 
ment. Let us assume that the system plus the environment starts from its i-th level. 
A straight path (e.g., green path in Fig. [TJ) still yields Plz because of (JT]). But now 
there are many other paths to excite the system, which in general may interfere 
with each other (e.g., blue and red paths in Fig. [TJ) . Those paths, however, can only 
add to the probability. Such property also exists if the environment is initially in 
a mixed state. Therefore, an initially excited (incoherent) environment will always 
increase the LZ transition probability. An important question now is to what extent 
the final ground state probability P g f will be suppressed as to grows. If P g f— >0 
as to— >oo, then AQC will be impossible in the presence of such an environment. 
Fortunately, as we shall show below, this is not the case. 

Phase coherence plays a role when different paths that start form a common 
level rejoin in another level and therefore interfere with each other. The phase that 
the wave function picks along each path is the time integral of the eigenenergy 
along that path. For very fast evolutions, all the paths come with the same phase 
(i.e., ~ 0). At slower evolutions, different paths will have different phases. Because 
of natural randomness in the parameters, the relative phases of these paths will 
be randomized and therefore will eliminate the effect of coherence. Thus, phase 
coherence is expected to have negligible effect at long times and large number of 
paths (i.e., large to). This is justified by our numerical simulations already at to — 4 
(see Fig. El). 

In practice, the environmental spins also are subject to thermal fluctuations, 
which affect the system via fluctuating local fields Bj. This effect can be modeled 
by coupling the spins to an additional (e.g. bosonic) environment which can cause 
additional dephasing between the energy levels making our incoherent approach 
even more relevant. Such an environment does not change LZ probabilities Pij in 
the fast passage limit (pij«l)P^ which is always the case when m— >oo. Here, we 
also neglect any relaxation effect as it can only reduce the excitation probability, 
i.e., enhance success probability in AQC. Hence, our calculation provides a lower 
bound for the AQC success probability. 



I 



I 
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Let otij and /3y be the probabilities for the upper and lower energy states just 
before traversing the ij-th anticrossing, respectively (see Fig. [IJ. We can write 

a(i-i) 3 = p lJ a lj + (1 - Pij) ^j, ^ 

Pi{j+X) = PijPij + (1 -Pij)&ij> 

(?) (i) (i) (i) 

with boundary conditions: 

<*Mj=P£', a 0j=Pff, Pu=P£', and A(M+i)=-P e / • Here 
Pfjp and Pj;p (Pj$ and P^)) are the initial (final) probabilities of finding the 
total system in |0, j) and respectively. This has been studied by switching to 

continuous variables by Ref. HH Here we instead solve the discrete relation while 
keeping the constraint ([11]) . From (JX3J) we get 

ai -i=P^^+L^ ai , 

Pj +1 =pVr,W+uWpj, (13) 
with vectors on=[om, ...,a iM ] T , Pj=\fiij, —,PMj] T , A*j — (1 - Pij)X^iZ\PiU 

Vi =(! - Pij)UZi+i py> and 

{p„ fc = j 

(1 - pa.)(l - Pi j) nCfc+i Pi* i > k ' 
fc > j 

!Pij i = k 

(l-Pk 3 )(l-p l3 )Ui=l + iPi 3 k> l ■ 
i > k 

We define the vectors P g i and P e i (P g f and P e /) to represent the initial (final) 
ground and excited state distribution, respectively. Here, we only consider the case 
where the central system is initialized in its ground state, i.e., P e i = 0, for which 
we have 

P ef = AP gi , with A = Z7Wt/( M - 1 )...[/W. (14) 

The final excited state probability P e f is found by summing over all elements of 
P e f, which is equivalent to tracing over the environmental degrees of freedom. 

Naturally, it is very difficult to find an analytical solution for the general case. 
Here, we consider a special case where all the gaps have the same size, gij=g m / ^/M, 
hence ptj = p = P^z ■ Such a situation occurs when \Jj\ = \Bj\. We find A = U M , 
where 

' p i=j 
Ua ■ = { p^-^l-p) 2 K] . (15) 
i > j 



Although U is not diagonalizable, we can still compute the M-th power using its 
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Jordan decomposition. Using induction, it can be shown that A=TJ M T~ 1 , where 



= 



Ti 



i j 



P M - j+i (^)i<j 

i > j ' 

1 j=i=i 
(-P) j - i (i-p)- 2 ^( j r_l)j>i>2 

else 

1 j=i=l 
j>i>2 . 

else 



After some manipulation, the final excitation probability becomes 

M 



Pef 
Vl = p 



M p 

3-2 



^>1 = P 



t=o ^ * 



2(t- 



■■CO*?' 



,t+s 



71/ 



(16) 



(17) 



where q = (1 —p)/p. Indeed, is the excitation probability when the environment 
is initially in the j-th state (Pg{ = 1)- Notice that for Plz > e _M , p » 1 and hence 
0<g<l. 




Fig. 2. vm as a function of Plz (= e _7r£ 'm'' 2 ''). The thin (black) lines are obtained using 11171 
with m = 0, 1, 7, ordered from bottom to top. The thick dashed (blue) and solid (red) lines are 
exact and approximate (M— >oo) solutions using (1 2 1 P and (122 ft . respectively. Inset: Vj vs Plz for 
m = 7 and, from bottom to top, j=l, 16,32,48,64,80,96,112, and 128. 
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The inset of Fig. [2] shows Vj vs Plz for different values of j, when m=7. The 
initial rapid rise happens when p (= P^ M ) changes from to ~1, thus becomes 
extremely sharp at large m. Notice that Vj monotonically increases as j grows. 
This reflects the fact that higher energy levels have more paths available to them 
for excitation. Maximum excitation happens when the system starts in the highest 
energy level (i.e., P^^ =1 for which P e f=VM), hence vm gives an upper bound for 
the excitation probability (P e f < %). Figure [2] displays vm as a function of Plz- 
As m grows, P e f deviates from Plz (m—0 case). The curves, however, saturate at 
m > 6. This is important because it shows that there is always a nonzero ground 
state probability even when m— >oo. All curves asymptotically join the m=0 curve 
as Plz— *-l, indicating that in the fast passage limit the environment does not affect 
the excitation probability. 

To find an analytical form for vm in the large M limit, we use (|17j) and write 
vm=Sq+Si, where 



M-2 

-,2(4+1) 



^^—E^ +i) ( M r 2 )(m 

§V(V)§<> (18) 

Using (-) = + 2(f- 2 ) + (It), we have 

^i = i[l+P 2M - 2 {2-p-*)W M ]-R M , 



2 

M 



m-2 , M _ 2> a 



t=0 

2A/-2 (M-2\ (M-2 



t=o 



2tH 



Now, using 



J=0 

which holds for x>0 and any closed contour C around the origin, one can show 
that So<pq<gil and i?M<(p?) 2 €l- Choosing an appropriate contour, (fT9|) yields 

W M = — (l + 2q cos 6> + g 2 ) ~ 2 d». (20) 

2 71 " J-7T 

AsM->oo, (1 + 2qcosO) M -► Pj^ cos6 . Hence 



if l+ r^. 

A-f^oo 2 \ ,/_„. 277 

For Plz <C 1, we can further approximate the integrand as a gaussian to get 

Voo » 1/2 + erf (vr V- In P L z)/4 ^-tt In P LZ - (22) 
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which gives Voo«l/2 in the adiabatic limit (Plz— >0). Equations (|2Tj) and ([22)) are 
plotted in Fig. [2j Equation (|22|) provides a good approximation for Vqo for a wide 
range of Plz- This has a very important implication for AQC. It shows that there is 
always a nonzero probability of success bounded from below by 1— Woo- In the adi- 
abatic regime (Plz^I), an environment of the sort modeled here can only suppress 
the probability of success by a factor close to 1/2. 

5. Numerical simulation 

We have also performed fully coherent simulations, using the Liouville equation 

p{t) = -i[H{t),p{t)], (23) 

for the density matrix p of the system plus the environment. For a case with 4 
environmental spins, plotted in Fig. [3j the coherent simulation shows only small 
fluctuations around the incoherent calculation, and also is very close to the equal 
gap calculation using (fT6|) . The figure also confirms that provides an upper 
bound for the excitation probability, as we expected. We have repeated this simu- 
lation for many other instances and obtained similar behavior. 

1 

0.8 
0.6 

ef 

0.4 



0.2 





0.2 0.4 0.6 0.8 1 

Fig. 3. Final excitation probability for completely coherent evolution and the approximate so- 
lutions. We assumed equally populated initial environmental states (P^'=l/M). For the "inco- 
herent", "equal gap", and ""Woo" curves we used Q14t . and l|17J l respectively. Parameters are: 
Sm =0.002, [B?,BJ]=[0.08,-0.04], [-0.07,0.10], [0.15,-0.11], [-0.23,0.06], Bj=0, and Jj=0.148 
for all m=4 environmental spins. 

6. Conclusions 

We have shown that a T—0 noninteracting spin environment does not affect the 
transition probability in adiabatic quantum computation in the small gap regime. 
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An excited (equilibrium or nonequilibrium) environment, on the other hand, will 
increase the excitation probability. We found that phase coherence does not play 
an important role in the excitation probability. Using an incoherent model with 
equal intermediate gaps, we identified a nonzero lower bound for the ground state 
probability. In the adiabatic limit, such a lower bound is ~l/2. AQC therefore is 
possible in the presence of such an environment. 
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Appendix A. Adiabatic Grover search problem 

In this appendix, we systematically derive the effective two-state Hamiltonian for 
the case of adiabatic Grover search (AGS) algorithm. In AGS, as defined by Roland 
and Cerf 1^01, the initial and final Hamiltonians are 

§ = i-|+)(+|, ^ = i-HR (A.i) 

where I is the unity matrix, |a) is the marked state, and 

with |Z) representing states in the computation basis. Let us define a new state 



which is orthogonal to |a). It is easy to see that |+) = (V7V — l\a) + \a))/y/N. The 
Hamiltonians (|A.1[) can therefore be written completely in the subspace made of 
|m) and \m): 



1 ( 1 \/N-l\ „ ._ _ fl 



Hi/E = I ~N\^/N=T N-l )' H ^ E = I -{00)- (A - 4) 
Using the Pauli matrices t x - z in the new 2-state subspace, 

H s = (l- s)Hi + sH f = - X -(lr z + ~g m r% (A.5) 

where 

. E+(N-l)e . ^EEl tw \ a c\ 

g = jj , 9m = — ^ — (E~e), (A.6) 

and e = E(2s — 1). We have thrown away trivial terms proportional to I. In the 
limit of large N, Hamiltonian (|A.5P reduces to ^ with g m — E/y/N. 
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Let us now introduce interaction between system qubits and environmental spins 
via the most general Hamiltonian 

n m 

= - E EEW 

a,/3=x,y,z i=l j=l 

Here &i and aj denote Pauli matrices that act on the i-th qubit and j-th environ- 
mental spin, respectively, and m is the number of the spins. In the subspace made 
of | a) and \a), one can write the reduced interaction Hamiltonian by substituting. 



^/N~T 

N - 1 V VN- 1 N -2 



1 / -is/W^l 



N - 1 \WN - 1 N-2 

1 

,0-l/(7V-l), 

Or in terms of Pauli matrices r 



(A.7) 



2{N - 1) 2(N-1) 
In the large N limit, one obtains 



N-2 N 



which yields 



where 



(A.9) 

V E E^ ( A -!0) 



2 

0=x,y,z j=l 



4 = E( J if- J if- J ff)- (A' 11 ) 

Since the environmental spins are uncoupled, we can arbitrarily define the 2-axis 
for each spin in such a way to obtain only longitudinal coupling: 



-ffint 

where 



in 

= -^E J ^< ( A - 12 ) 



4 = ^/;) 2 + g/]0 2 + g/;) 2 - (A.i3) 

This result is valid more generally than for the spin environment studied here. In 
other words, at least for the adiabatic Grover search problem, any environment that 
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couples to the individual qubits in a most general way would result in a longitudinal 
coupling of the effective two-state system (made of the lowest two energy states) 
to the environment. This can be clearly seen from Eqs. (|A.8|) and (|A.9|) that are 
independent of the type of environment. Notice that in (|A.8|) the terms involving 
t x ' v always appear with a coefficient of 0(1/ ^/N) and therefore will vanish in 
the limit N — ► oo. It physically makes sense because such off-diagonal terms in 
the interaction Hamiltonian would cause relaxation between the two levels with 
a rate of 0(1/ N). This means solving the search problem just based on classical 
relaxation to the ground state (i.e., classical annealing) will take a time of O(N) 
which is the complexity of the classical search. One therefore would expect that the 
same property also persist for problems that are more general than the adiabatic 
Grover search; otherwise, e.g., spin glass problems could be solved efficiently by 
just letting the system relax to the ground state. 
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